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ABSTRACT 

We study the dispersal of a gaseous disc surrounding a central high-mass stellar core 
once this circumstellar disc becomes fully ionized. If the stellar and surrounding EUV 
and X-ray radiations are so strong as to rapidly heat up and ionize the entire circum- 
stellar disc as further facilitated by disc magnetohydrodynamic (MHD) turbulence, a 
shock can be driven to travel outward in the fully ionized disc, behind which the disc 
expands and thins. For an extremely massive and powerful stellar core, the ionized gas 
pressure overwhelms the centrifugal and gravitational forces in the disc. In this limit, 
we construct self-similar shock solutions for such an expansion and depletion phase. 
As a significant amount of circumstellar gas being removed, the relic disc becomes 
vulnerable to strong stellar winds and fragments into clumps. We speculate that disc 
disappearance happens rapidly, perhaps on a timescale of ~ 10 3 — 10 4 yr once the 
disc becomes entirely ionized sometime after the onset of thermal nuclear burning in 
a high-mass stellar core. 



Key words: circumstellar matter 
early type - stars: winds, outflows 



hydrodynamics - plasmas - shock waves - stars: 



1 INTRODUCTION 

There are ample evidence that young stars are often ac- 
companied by circumstellar discs, as part of star formation 
processes. Discs are most frequently detected for low-mass 
T Tauri protostars. The probability of detecting a circum- 
stellar disc decreases with an increasing stellar mass. Natta 
et al. (2000) concluded that a lack of discs around Herbig 
Be stars may result from their faster evolution. Fuente et 
al. (2003) revealed the first evidence of dusty discs around 
Herbig Be stars, yet the disc-to-star mass ratio Md/M* is at 
least an order of magnitude lower than those of Herbig Ac 
and T Tauri stars, suggesting that discs may evolve much 
more rapidly around very massive stars. Suppose that high- 
mass stars form similarly as low- or intermediate-mass stars 
in the sense of involving collapse (e.g., Lou 1996; Lou & 
Shen 2004; Bian & Lou 2005; Yu & Lou 2005; Yu, Lou, 
Bian & Wu 2006; Lou & Gao 2006) and accretion (e.g., Mc- 
Kee & Tan 2002; Blum et al. 2004), then a circumstellar disc 
should form during the early evolution phase. Observational 
evidence for accretion and massive discs during high-mass 
star forming stages were indeed reported (e.g., Sandell et al. 
2003; Pestalozzi et al. 2004). While massive stars of late O or 
early B types may have disc signatures, more massive stars 



lack clear evidence for discs (e.g., Blum et al. 2004). Unless a 
massive star forms through other processes without involv- 
ing a circumstellar disc, the key issue is to understand the 
mechanism for a fast and efficient dispersal (primarily hy- 
drogen gas) of circumsteilar discs around the most massive 
stars. 

Circumstellar discs are important because they are sites 
where a companion star or planets may form. For a massive 
star, a considerable mass fraction of its dense circumstel- 
lar disc appears to be rapidly removed. The currently esti- 
mated lifetime of circumstellar discs for low-mass stars is of 
the order of several million years (e.g., Strom et al. 1993; 
Bricefio et al. 2001). For intermediate- and high-mass stars, 
the timescale for the existence of circumstellar discs can be 
much shorter because of faster evolution and stronger radia- 
tive power output from central stars (e.g., Natta et al. 2000). 
Theoretically, the dispersal or destruction of a circumstellar 
disc may be modelled through viscous disc accretion onto a 
central star, photoevaporation by either stellar or external 
UV radiations, interactions with other stars or consumption 
by planet formation, or any combinations of these processes 
(e.g., Hollenbach et al. 2000). The photoevaporation and ac- 
cretion models have been widely applied to disc dispersals of 
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low-mass T Tauri stars (e.g., Johnstone et al. 1998; Storzer 
& Hollenbach 1999; Clarke et al. 2001). 

We here advance a plausible physical scenario for the 
dispersal of circumstellar discs around most massive stars. 
By forming a massive protostellar core, the accretion gradu- 
ally ceases and a massive (perhaps > 0.1 — 0.3M Q ) centrifu- 
gally supported H I disc 1 forms around the core. Near the 
main-sequence when the stellar core ignites thermal nuclear 
reactions to produce a profuse amount of extreme-ultraviolet 
(EUV) photons, an ionization front (IF) travels rapidly out- 
ward to ionize the entire H I disc to a high temperature 
> 10 4 K (e.g., Osterbrock 1989). Meanwhile, X -ray emis- 
sions associated with protostellar and disc magnetic activi- 
ties also ionize gas to ~ 10 6 — 10 7 K within a few thousand 
AUs (e.g., Glassgold et al. 2000). For a high-mass star of 
M* ~ 10M Q , the stellar luminosity alone may not fully ion- 
ize a very dense circumstellar disc. However, a circumstellar 
disc may be fully ionized if (a) the central stellar mass M* is 
much greater than lOM© or the Lyman flux is considerably 
higher; (b) the host star is born in a cluster, where its cir- 
cumstellar disc is exposed to intense external X-ray and UV 
radiations; (c) a circumstellar disc has been diluted signif- 
icantly via other processes such as photoevaporation (e.g., 
Hollenbach et al. 1994); (d) magnetorotational instabilities 
(e.g., Balbus & Hawley 1998; Balbus 2003) in the weak mag- 
netic field regime sustain a level of MHD turbulence to inces- 
santly mix ionized and neutral components and to effectively 
enhance ionization. To estimate, we assume a diluted neutral 
Hi disc 2 with a surface mass density E(r) = Ed(r / R)^ 1 ^ 2 . 
For M d ~ O.O2M , we derive E d = 3M d /(47r.R 2 ) ~ 
0.042 g cm~ 2 such that the H I column density normal to 
the disc is N H (r) ~ 2.5 x W 22 (r/R)- 1/2 cm" 2 . The Ly- 
man flux incident on the disc will penetrate a column 
N e (r) < 1.8 x 10 22 $^ 2 (r/i?)- 1/2 cm" 2 (see Appendix A 
of Hollenbach et al. 1994). Thus an H I disc may be fully 
ionized 3 for $49 > 2, which is readily satisfied for early 
O stars (e.g., Panagia 1973). For more massive stars with 
$49 ~ 10 — 100, a denser circumstellar disc of mass Md 
even higher than ~ 0.2M Q might be ionized completely as 
facilitated by X-ray fluxes and disc MHD turbulence. If sub- 
merged in an external X-ray and EUV radiation field from 
nearby massive stars within a local star cluster (which is 
the usual case), a complete ionization of a circumstellar disc 
becomes more likely. A hot H II disc inevitably expands as a 
result of the overwhelming pressure gradient. Typically, the 
gas density is higher in the inner part than in the outer part 
of a disc, so does the pressure gradient. Hence, the inner disc 



1 Discs around massive stars may be even more massive initially 
(see Sandell et al. 2003). 

2 Disc radii R of high-mass stars are expected to be larger than 
those of T Tauri or Herbig Ae/Be stars, which extend from tens 
to hundreds of AUs (e.g., McCaughrean 1997; Meyer & Bcckwith 
2000). We take R ~ 10 3 AUs (e.g., Pestalozzi ct al. 2004). 

3 Around massive stars, dust grains are effectively evacuated to 
larger radii (~ 0.1 pc) due to sublimation, stellar winds and ra- 
diation pressure (e.g., Chini et al. 1987; Churchwcll 2002). When 
dust attenuation is included, the required relic disc mass will be 
smaller. In the approximation of Franco ct al. (1990), the inci- 
dent Lyman flux $49 is reduced by a factor e~ Td where r d is the 
optical depth at the boundary of HII regions. For -r d > 1, the disc 
mass that can be fully ionized is reduced by a factor less than 0.6. 



expands faster than the outer disc does and consequently a 
shock develops in the disc to travel radially outward. The 
disc portion behind the shock will continue this expansion 
to reduce the disc surface mass density. 

The main thrust of this Letter is to construct a gas 
dynamic model for a such disc expansion phase after a cir- 
cumstellar disc has been ionized and to estimate disc disap- 
pearance timescales using available data. While our results 
are preferably applicable to high-mass stars ( > 10M©), 
they are also relevant for discs around low- or intermediate- 
mass stars in binaries where the disc of the secondary may 
be fully ionized by the X-ray and EUV field of the more 
massive primary. 



2 DISC MODEL FORMULATION 

We start with a specific case of a neutral H I disc carrying a 
surface mass density E oc r _1 , and then extend the results to 
different power-law indices of E profiles (Shen & Lou 2004a; 
Shen, Liu & Lou 2005; Lou & Bai 2006; Wu & Lou 2006) in 
the limit when the thermal pressure dominates other forces 
in the shock evolution of an H II disc. The stellar gravity 
has weaker effects over distances larger than f ~ GM*/(2a 2 ) 
where a is the sound speed of the H II gas. For M» ~ 10M Q 
and a ~ 10 km s" 1 , we have f ~ 40 AU and may neglect 
the stellar gravity on the main portion of the disc with a 
typical radius R ~ 1000 AU. 

We presume a thin neutral Hi circumstellar disc around 
a protostar in a rotational equilibrium initially with a sin- 
gular isothermal disc (SID) profile, namely 



E(r) = 



a ( 2 ,(l + D 2 
2?rGr 



M{r)= " y ^ ' , j(r)=a Dr, 

(1) 

where in cylindrical coordinates (r,8,z), E(r), M(r) and 
j(r) are the surface mass density, the enclosed mass inside 
radius r and the 2— component specific angular momentum, 
respectively; G is the gravitational constant and D is the 
rotational Mach number defined by D = Vo/ao with Vo and 
ao being the rotation speed 4 and the isothermal sound speed 
of the H I gas. 

At a certain initial time t = 0, thermal nuclear reactions 
start in the stellar core and the resulting radiation rapidly 
ionizes the circumstellar disc in a short time to reach an 
isothermal state with a faster sound speed a. To keep the 
initial SID profile unchanged would require 



S(r, 0+) = eV (1 + g2) , j(r, 0+) = eaDr 



M(r,0 + ) = 



27rGr 
e 2 q 2 (l + D 2 )r 
G 



(2) 



where e = ao/a < 1 is the square root of the temperature 
ratio. Now the pressure and centrifugal forces overwhelm 
the gravitational force and a circumstellar disc expands in- 
evitably. As the pressure gradient is steeper in the inner disc, 
a shock will emerge to travel outward. This is a close analog 



4 We take on a flat rotation curve instead of the usual Kcplcrian 
rotation curve. As can be seen later, the centrifugal force can be 
omitted in the limit when the pressure dominates and thus the 
choice of different rotation curves does not really matter. 
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of the so-called 'champagne flows' in H II clouds surround- 
ing OB stars (Tenorio-Tagle 1979; Franco et al. 1990). Since 
the latter may evolve towards a self-similar phase (Shu et 
al. 2002), it is suggestive of a similar description for a disc 
expansion with a shock. Here, we focus on axisymmetric 
similarity shocks relevant to the disc dispersal scenario. 

Under the axisymmetry, we use the following nonlinear 
equations to describe disc dynamics in cylindrical geometry. 



Table 1. Parameters for numerical solutions 
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where u is the radial speed, <!> is the gravitational potential, 
II is the two-dimensional thermal gas pressure, M(r, t) = 
J Q r E(r' ,i)2nr' ' dr' is the enclosed mass within radius r. The 
Poisson integral is 



$(r,t) = 



-GE(r',t)r'dr'd0 



( r /2 +r 2 _ 2r'rcos8) 1 / 2 
We introduce the following similarity transformation 
aa(x) 



r 
at 



E(r,t) 



2ixGt 



M(r,t) = —m{x) 

Gr 



w(r, i) = av(x) 
j(r,t) = Pa 2 tmix) 



(6) 



(7) 



in the two-dimensional isothermal self-gravitating nonlin- 
ear ideal gas equations — JBJ, where a(a:), v(x) and m(x) 
are the reduced surface mass density, radial flow speed and 
enclosed mass, all being dimensionless functions of x. Pa- 
rameter (3 = D/[e(l + D 2 )] is determined from the equilib- 
rium. The last one in transformation Q holds as the ra- 
tio j(r)/M(r) remains constant since t = by the angular 
momentum conservation (Li & Shu 1997; Saigo & Hanawa 
1998). 

In the so-called monopole approximation, 5 the similar- 
ity ordinary differential equations (ODEs) become 



[{x - v) 2 - 1] dv _ a(x -v)[l- (3 2 aix ■ 



(x — v) dx 



\{x - v) 2 - 1] da _ a[l - (3 2 a(x - v)] - (x - v) 



(x — v)a dx 



(8) 



(9) 



which, together with m = x(x — v)a, are to be solved for 
specified 'asymptotic conditions'. 

The leading asymptotic solutions at x — - > +oo for either 
t — *• + or r —> oo are 



v -» V + 



1 - A + f3 2 A 2 



A 

a — > — 

x 



m^Ax, (10) 



where V and A are two constants and higher-order terms 

5 A multipole treatment can be found elsewhere (Li & Shu 1997). 
This monopole approximation does not affect our results much 
because when £<1, the self-gravity becomes negligible (Shu et 
al. 2002; Shen & Lou 2004b). 
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The three rows without values for D and e arc unphysical cases. 

can also be readily derived. With an initial SID, boundary 
conditions for x —* +oo are V = and A — e 2 (l + D 2 ). 

There exist two classes of asymptotic solutions at x — * 
for either r — > or t ^ oo; we focus on the class of regular 
solutions with leading behaviours of 



v(x) — > x/2 



a(x) — > B , m(x) — » Bx 2 /2 (11) 



as x — > 0, where B is an integration constant. 

The two isothermal shock jump conditions are 

\2 



(v d - X s )(v u 



= i 



a, 



= ivd - x s ) a d , (12) 



where subscript d (u) denotes the downstream (upstream) 
of a shock and x s is the 'shock location' in the similarity 
variable x with a shock speed ax s (Shen & Lou 2004b; Bian 
& Lou 2005; Yu et al. 2006). 

To solve ODEs JSJ and JUJ numerically, we specify /3 
leading to a relation for D and e by f3 = D/[e(l + D 2 )]. A 
shock solution is constructed as follows: (i) Integrate ODEs 
JSJ and from the origin using solution 11 U for various 
B; (ii) Impose shock condition 11211 at each integration step 
to cross the sonic critical line x — v = 1 and continue for 
larger x; (iii) Pick out the one that matches condition 1101 as 
x — > oo with V = 0. By matching a proper A value, this last 
step leads to another relation of D and e by A — e 2 (l + D 2 ). 
Both D and e are then determined accordingly. 

Remarkably, it is found that in the limit B — ► for a 
disc heated to high temperatures, there exists an invariant 
shock solution not affected by /3. The A value is uniquely 
determined by B with A — ■> as B — > 0, while the ratio A/B 
approaches a constant ~ 1.58 (Table 1). Physical solutions 
require A/3 2 < 1. This bears a strong resemblance to the 
spherical similarity solutions of 'champagne flows' (Shu et 
al. 2002). This invariant form describes a situation when 
the thermal pressure overwhelms in the disc expansion with 
negligible gravity and centrifugal forces. The limit of x 3 = 
2.00 gives a shock speed twice the ionized sound speed a. 
The downstream (post-shock) gas becomes more uniformly 
distributed and a linear expansion emerges from small x. 

Following the above procedure by ignoring gravity and 
centrifugal forces, one can extend these results to a power- 
law surface mass density E(r) oc r~ n where < n < 2 (Shen 
et al. 2005). Similar to the n = 1 case, the downstream 
behind the shock approaches a uniform density yet with the 
enclosed mass oc t~ n and the flow speed becomes linear in 
radii, viz., v — > nx/2 as x — > 0. 
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Figure 1. A shock solution for v(x), a(x) and mix) versus x. The 
dash-dotted line in the top panel is the sonic critical line x — v = 1. 
The reduced radial speed v(x) and density a(x) encounter discon- 
tinuities at shock location x s = 2.00, while the reduced enclosed 
mass m(x) is continuous by the mass conservation and 1121 . 
Downstream (post-shock) and upstream (pre-shock) portions arc 
shown in heavy and light solid curves; the dotted curves denote 
the leading-order terms of asymptotic solutions 1111 . 



3 CIRCUMSTELLAR DISCS 

We use our self-similar invariant shock solutions to esti- 
mate a circumstellar disc dispersal timescale. To elucidate 
the concept, we introduce the following fiducial solution us- 
ing the SID model. For a centrifugally supported neutral 
equilibrium SID, we may choose D = 2 for a supersonic 
rotation against self-gravity 6 . The Lyman ionization may 
lead to a ao/a ratio of e = 0.02 (e.g., ao ~ 0.2 km s _1 and 
a ~ 10 km s _1 ). We then have A — 2 x 10~ 3 which gives a 
B = 1.26 x 10~ 3 and a (3 = 20. Yet this solution remains 
essentially the same as that of (3 — and B = 1.26 x 10 -3 . 
That is, there exists an invariant shock solution (Fig. 1) not 
affected by j3 in the limit of B — > with a shock at x s = 2.00. 

Since a similarity solution does not carry characteris- 
tic timescales or lengths, one needs a starting point to ap- 
ply the solution in proper spatial and temporal domains. 
For example, at a particular radius R ~ 1000 AU, i.e., the 
ionized disc radius after some time of expansion, the corre- 
sponding similarity variable is x c — R/(at c ) where t c is an 
input timescale that a self-similar evolution has lapsed and 
marks the starting epoch for the subsequent evolution. Here, 
t — t c is a calibration point. We focus on R and presume 
that for t > t c , all materials leaving R would be regarded 
as lost from a circumstellar disc. The key issue of concern is 
the evolution of the enclosed disc mass M(R, t) for t > t c 
where t c is a starting time when the mass and radius of a 
disc are specified. The information may be extracted from 
the bottom panel of Fig. 1, i.e., the downstream portion of 
m(x) in heavy solid curve. For a small x < 0.5, the reduced 
mass is m(x) = Bx 2 /2 and the evolution of M(R, t) may 



6 D should be > 1, but may not be too large. Our basic results 
derived here are not very sensitive to D variations. 



be estimated by 

M{R, t) aBR 2 /(2Gt) . (13) 

As a self-similar evolution starts from t ~ 0, it takes a 
timescale of t c = R/(2a) ~ 250 yr for a shock to reach R, 
where t c is the calibration time corresponding to a similarity 
coordinate x c = 2. At this epoch, the disc mass behind a 
shock is M(R,t c ) = a 3 t c ml2)/G ~ 0.2M Q where m(2) = 
3.4 x 10 -3 by the heavy solid curve of m(x) in Fig. 1. From 
now on, the shock continues to travel into distant regions and 
eventually merges into ambient Hll clouds; meanwhile, the 
subsequent self-similar evolution of the disc is well described 
by the downstream (post-shock) portion. We may then use 
interpolated values of m(x) with x < 2 and M{R, t) = 
a 3 trn(x)/G to estimate M(R, t) for t > t c . For a sufficiently 
small x or large t, we may directly use equation 11311 that 
shows an enclosed disc mass proportional to t . At later 
times of 5t c (x = 0.4), 10t c {x = 0.2), 20t c (x = 0.1), 100t c 
(x = 0.02), corresponding to ~ 1250 yr, ~ 2500 yr, ~ 5000 
yr, ~ 25000 yr, the enclosed disc mass becomes ~ 0.03Mq, 
~ O.O15M , ~ O.OO75M , ~ 0.0015M Q , respectively. 

Diagnostically, it is not easy to set a criterion for 
the 'disc disappearance' by decreasing enclosed disc mass 
M(R,t). We propose tentatively that a 'disc disappearance' 
corresponds to a disc mass dropped between one tenth and 
one hundredth of the mass ntt = t c . This is partly due to the 
fact that a disc depletion leads to fainter Ha emissions from 
ionized gas and infrared/submillimeter emissions from dusts, 
and partly due to the possibility that a rarefied disc be- 
comes unstable by strong stellar wind shears and breaks into 
clumps Thus, a timescale of disc disappearance is estimated 
as tdis ~ 10t c — 100t c by expression 113H . For t c ~ 250 yr, the 
'disappearance timescale' of tdis ~ 2.5 x 10 3 — 2.5 x 10 4 yr is 
very much shorter than the typical timescale of several Myr 
for low-mass T Tauri stars. 



4 DISCUSSIONS 

While the self-similar process described here is highly ideal- 
ized, it may catch several gross features of a disc dispersal 
process not fully explored so far. First, we assume a fast 
traveling IF together with environmental X-ray and EUV 
radiation field that heat up the entire disc to an isother- 
mal state in a short time. Secondly, rotation and gravity are 
ignored in the limit of an overwhelming thermal pressure. 
Our model of invariant shock is more applicable to discs 
around very massive stars or to circumstellar discs exposed 
to intense X-ray and EUV radiation from nearby external 
sources (e.g., evolved massive stars or supernova explosions; 
Chevalier 2000). Thirdly, we ignore the disc thickness. In re- 
ality, such a disc should expand both radially and vertically 
besides rotation, with photoevaporation being concurrent. 
A thorough analysis of these aspects deems worthwhile. 

We now summarize highlights of our model. First, a cir- 
cumstellar disc as a whole is expanding rather than fixed or 
shrinking shortly after the ignition of nuclear reactions in 
the stellar core. This expansion induced disc mass-loss rate 
is modelled as self-similar and differs from classic photoe- 
vaporation models. For very massive stars submerged in an 
intense X-ray and EUV radiation field, circumstellar discs 
become too short-lived to be observed after entering the 
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main-sequence. Secondly, the fast moving shock serves to 
disperse disc materials and the rarefied disc may become 
vulnerable to fragmentations by stellar wind shears. These 
processes accelerate the disc destruction. Thirdly, this self- 
similar expansion tends to establish a more uniform density 
distribution even though the initial density profile is non- 
uniform. Fourthly, during this self-similar expansion phase, a 
considerable mass fraction of the disc is expelled or displaced 
to distant places; these processes might explain the presence 
of appreciable amounts of ionized gas at intermediate radii 
(i.e. 10 3 ~ 10 4 AU) in ultracompact H II regions (e.g., Wood 
& Churchwell 1989; Shu et al. 2002). Moreover, as time goes 
on, these dispersed disc materials eventually cool down and 
may become a reservoir for producing Edgeworth-Kuiper 
Belt-like objects around a central star. 

We now comment on several aspects of further model 
development. First, even for an isothermal model, the tem- 
peratures across a shock can be allowed to be different (Shen 
& Lou 2004b; Bian & Lou 2005; Yu et al. 2006). Secondly, 
the isothermal condition can be replaced by the more gen- 
eral polytropic approximation (Wang & Lou 2006; Lou & 
Gao 2006). Thirdly, it is of considerable interest to incorpo- 
rate the effect of a magnetic field (Shen et al. 2005; Yu & 
Lou 2005; Lou & Zou 2004, 2006; Lou & Wu 2005; Wu & 
Lou 2006; Lou & Bai 2006). 

Finally, we speculate that grossly spherical 'champagne 
flows' with shocks traveling in H II clouds (e.g., Franco et al. 
1990; Shu et al. 2002; Shen & Lou 2004b; Bian & Lou 2005) 
on much larger spatial scales encompassing a circumstellar 
disc can be driven by an intense central stellar radiation, 
stellar winds and outflows from the disc sustained by photo- 
evaporation. Triggered by shocks and various flow or ther- 
mal instabilities, fragmentations would occur in clouds; by 
cooling and coagulation with time, these clumps may evolve 
into comet-like objects. This might be the origin of the Oort 
cloud around our solar system and should be common in 
other protostellar systems as well. 

Not only for young stellar objects, circumstellar discs 
are also observed around main-sequence stars or even more 
evolved stars (e.g., Zuckerman 2001). These debris dusty 
discs are exposed to intense radiation from central objects 
during the late evolutionary stage. Disc expansions and 
shocks can be initiated in these discs and evolve in a self- 
similar manner. Even with considerable idealization, the 
analysis is technically challenging as rotation, central grav- 
ity and disc self- gravity need to be taken into account. 



ACKNOWLEDGMENTS 

This research was supported in part by the ASCI Cen- 
ter for Astrophysical Thermonuclear Flashes at the U. of 
Chicago under Department of Energy contract B341495, by 
the Special Funds for Major State Basic Science Research 
Projects of China, by the THCA, by the Collaborative Re- 
search Fund from the National Natural Science Founda- 
tion of China (NSFC) Young Outstanding Overseas Chi- 
nese Scholars (NSFC 10028306) at the National Astronom- 
ical Observatories, Chinese Academy of Sciences, by NSFC 
grants 10373009 and 10533020 at the Tsinghua Univ., and 
by the SRFDP 20050003088 and the Yangtze Endowment 



from the Ministry of Education at the Tsinghua Univ. Af- 
filiated institutions of Y-QL share this contribution. 



REFERENCES 

Balbus S. A., 2003, ARA&A, 41, 555 

Balbus S. A., Hawley J. F., 1998, Rev. Mod. Phys., 70, 1 

Bian F. Y., Lou Y.-Q., 2 005, MNRAS, 363, 1315 

Blum R. D., et al., 2004, |astro-ph/040919'o| 

Briceno C, et al., 2001, Science, 291, 93 

Chini R., Kriigel E., Wargau W., 1987, A&A, 181, 378 

Chevalier R. A., 2000, ApJ, 538, L151 

Churchwell E., 2002, ARA&A, 40, 27 

Clarke C. J., Gendrin A., Sotomayor M., 2001, MNRAS, 328, 485 
Franco J., Tenorio-Tagle C, Bodenheimer P., 1990, ApJ, 349, 126 
Fuente A., et al., 2003, ApJ, 598, L39 

Glassgold A. E., et al., 2000, Protostars and Planets IV, eds. V. 

Mannings, A. P. Boss, S. S. Russell (Tucson: Univ. of Arizona 

Press), pp 429-455 
Hollenbach D., Johnstone D., Lizano S., Shu F. H., 1994, ApJ, 

428, 654 

Hollenbach D., et al., 2000, Protostars and Planets IV, eds. V. 

Mannings, A. P. Boss, S. S. Russell (Tucson: Univ. of Arizona 

Press), pp 401-428 
Johnstone D., Hollenbach D., Bally J., 1998, ApJ, 499, 758 
Li Z.-Y., Shu F. H., 1997, ApJ, 475, 237 
Lou Y.-Q., 1996, MNRAS, 279, L67 
Lou Y.-Q., Shen Y., 2004, MNRAS, 348, 717 
Lou Y.-Q., Bai X. N., 2006, MNRAS, submitted 
Lou Y.-Q., Gao Y., 2006, MNRAS, submitted 
Lou Y.-Q., Wu Y., 2005, MNRAS, 364, 475 
Lou Y.-Q., Zou Y., 2004, MNRAS, 350, 1220 
Lou Y.-Q., Zou Y., 2006, MNRAS, 366, 1037 
McCaughrean M. J., 1997, ASP Conf. Ser. ed. David Soderblom, 

119, 53 

McKee C. F., Tan J. C, 2002, Nature, 416, 59 
Meyer M. R., Beckwith S. V. W., 2000, Lecture Notes in Physics, 
548, 341 

Natta A., et al., 2000, Protostars and Planets IV, eds. V. Man- 
nings, A. P. Boss, S. S. Russell (Tucson: Univ. of Arizona 
Press), pp 559-588 
Osterbrock D. E., 1989, Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei (Mill Valley: University Science Books) 
Panagia, N. 1973, AJ, 78, 929 
Pestalozzi, M. R., et al., 2004, ApJ, 603, L113 
Saigo K., Hanawa T., 1998, ApJ, 493, 342 
Sandell G., Wright M., Forster J. R., 2003, A pJ, 590, L45 
Shen Y., Lou Y.-Q., 2004a, MNRAS, 353, 249 jastro-ph/0405444) 
Shen Y., Lou Y.-Q., 2004b, ApJ, 611, L117 jastro-ph0407328t 
Shen Y., Liu X., Lou Y.-Q., 2005, MNRAS, 356, 1333 
Shu F. H, et al., 2002, ApJ, 580, 969 
Storzer H., Hollenbach D., 1999, ApJ, 515, 669 
Strom S. E., et al., 1993, Protostars and Planets III, pp 837-866 
Tenorio-Tagle G., 1979, A&A, 71, 59 
Wang W. G., Lou Y.-Q., 2006, MNRAS, submitted 
Wood D. O. S., Churchwell E., 1989, ApJS, 69, 831 

Wu Y., Lou Y.-Q., 2006, MNRAS, submitted 

Yu C, Lou Y.-Q., 2005, MNRAS, 364, 1168 I astro-ph0509242 I 
Yu C, Lou Y.-Q ., Bian F.-Y., Wu Y., 2006, MNRAS, in press 

lastro-ph/0604261 1 
Zuckerman B., 2001, ARA&A, 39, 549 



